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Abstract 

Thermal conduction across a magnetic field is 
strongly suppressed compared with conduction along 
the field. However, if a fiare is heated by a highly 
filamented beam directed along the field, then the array 
of heated cells in a cross-section of the fiare will result 
in both small spatial scales (with consequently large 
temperature gradients) and a large surface area for the 
heated volume, providing a geometrical enhancement 
of the total cross-field energy fiux. To investigate the 
importance of this filamentary geometry, we present a 
simple model of a single heated filament surrounded by 
an optically thin radiating shell, obtain an analytical 
expression for the stable equilibrium temperature profile 
within the shell, and use this to impose limits on the size 
of filament for which this model is appropriate. 

We find that this mechanism by itself is capable of trans- 
porting a power of the same order as a large flare, with 
a moderate range of filament sizes. The length scales are 
substantially smaller than can be resolved at present, al- 
though they should be regarded as underestimates. 
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1 Introduction 



seems to be a theoretical consensus both that prominences 



There is an observationa l consensus ( Engvold, 1976 ; Chi- 
uderi Drago et al., 1992 ) that solar prominences have a 
good deal of fine structure within them, on scales ranging 
from 100 km to 1500 km; similarly, each improvement in 
resolution finds more structure within fiares. Other stud- 



ies, such as Fletcher and Brown (1995), provide indirect 
support for the existence of fine filamentary structure in 
flares, by invoking the presence of such structure to ex- 
plain polarisation observations. At the same time, there 



condense due to a thermal instability (Field, 1965) in the 
plasma around them, and that cross-field thermal con- 
duction plays a significant role in determining the scale 
of their fine structure (Van der Linden, 1993). However, 
these studies (as reviewed by (Heyvaerts, 1974)) concen- 
trate on the dynamics of the instability which forms the 
threads, rather than on the final equilibrium of the result- 
ing structure. In this paper, we instead present a simple, 
quasi-static, analytical model of a filamentary structure 
in a prominence or fiare, which exhibits the small spatial 
scales necessary to produce a temperature gradient large 
enough to make cross-field thermal conduction significant, 
and which potentially gives cross-fleld conduction a dom- 
inant role in governing both the number and size of the 
threads in a prominence or flare, and the energy budget 
of the beam or current supplying the thread. 

As discussed in a preli minary version of this 
model ( Brown and Gray, 1994 ), we have available a large 
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number of free parameters, so the aim has to be not to 
make a single prediction, but rather to see which combi- 
nations of parameters allow a consistent model. 

The paper is organised as follows: in section ^ we de- 
scribe the model we have used; in section ^ we describe 
the analytical and numerical solutions to this model that 
we have obtained; in section 0we discuss the stability of 
these solutions; and in section]^ we briefly discuss the mi- 
nor changes to our results which are necessary in the case 
of a low-/? plasma. 

2 The model 

We suppose the existence of a heating mechanism which is 
spatially fragmented with fllaments aligned along the am- 
bient magnetic fleld, and so divides the heated plasma into 
a number of identical cylindrical cells, radius R, each of 
which is heated only within some coaxial cylindrical core 
(we have also considered the case of filamentation into 
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plane sheets, with similar results to those presented here). 
This heated core is surrounded by a shell which is heated 
solely by conduction from the core, and cooled solely by 
radiation. We find the temperature structure of this shell 
which allows an equilibrium to form by radiating away all 
the energy that is transported into it. We can ignore the 
effects of longitudinal conduction, and assume that the fil- 
amentation allows the transverse conducti on t o dominate 
- this is retrospectively justified in Sect. |3.4 Thus, the 
equilibrium structure of a particular 'slice' (of length /) 
through the cell will be governed only by the power input 
at a particular point along the cylinder. When we wish to 
find numerical values, below, we will take I = 10^ m as the 
cylinder length (this is the approximate stopping distance 
for a typical electron beam). 

We have been deliberately vague about the heating 
mechanism - we will find that the only relevant parameter 
is the total power per unit length conducted from within 
the core, so that this heating mechanism could equally 
well be a non-thermal beam or a dissipating current. The 
cost of this generality is that we can say nothing about 
the structure of the core, nor fix what fraction of the total 
power available is conducted out, rather than being car- 
ried longitudinally within the core. We expect, however, 
that the heat-loss from the core will be substantial. 

We parameterise distances from the centre of the cell 
in terms of the dimensionless parameter s — r/R, where 
s = 1 at the edge of the cell. We will take the shell to 
consist of material below some temperature Tp which we 
will take, below, to be around the peak in the Cox- Tucker 
radiative loss function. This shell starts a distance pR 
from the centre of the cell, so that p is defined by 



T(s = p)=Tp 



Conduction need not be the only method by which the 
core is cooled, but we will denote by Pp the power ac- 
tually conducted into the shell per unit length along the 
field. Finally, we take each cell to be independent, so that 
there is no energy flux between cells; thus the temperature 
gradient must vanish at s = 1. 

We emphasise that all the power input is confined within 
the radius p, and that we will not concern ourselves further 
with the details of this input nor the internal structure 
of the core. For this to be valid we must, for example, 
take the plasma in the shell to be optically thin so that 
it will not experience significant secondary heating by any 
radiative emission from the core. 

In the shell, there will be two processes competing: 
there will be a conductive heat flux entering from the 



core, and radiative emission from the heated plasma. 
For p < s < 1, the conductive heat flux per unit axial 
length is 

P{s) = Q{s)2ttsR, (2) 

where the heat flux density Q{s) is (neglecting K||dT/dr 
along the beam) 



ar 



(3) 



and the cross-field thermal conduction coefficient is 



(Rosenbluth and Kaufman, 1958) 



^22^1/2 



(4) 



with rt in m i? in Tesla, and 

Ki = 9.76 X 10"^^ W m^ T^ K"^/^ 

(note that the effective coefficient is likely to be rather 
larger than this - see the discussion in section |3.2| below). 
High in the solar atmosphere, magnetic pressure domi- 
nates, and we take the magnetic field to be approximately 
the same across the filament, so that the pressure-balance 
condition reduces to 

n{s)T{s) — a, (constant) (5) 

(but see Sect. ^) and in terms of kq = ki/B'^ we can write 

n^ = no{B)^ = Ko{B)a'T~'/\ (6) 



Given that the emitted power from a volume dV^ is 
(1) f{T)dV , for a given annulus at radius s, we will have 



P{s) - P{s + ds) - n{s)^f{T{s)) ■ 2ttsR ■ Rds = 0, 



which, on substitution of Eqn. 



reduces to 



— [sQ(s)]+7i2/(T)si? = 0. 
ds 



(7) 



The vanishing of the radial flux at the outer boundary of 
the slice provides one boundary condition (5(1) = 0, and 
another will be provided indirectly by the specified total 
power input to the shell at radius p. 

Before going on to determine the temperature profile ex- 
plicitly, it is worthwhile to point out that the forms of the 
transverse conduction coefficient (Eqn. (^)) and radiative 
loss function /(T) conspire to produce a temperature pro- 
file which is radically different from the profile produced 
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by longitudinal conduction (as in the usual treatment of 
flares and the transition region). Because kj_ depends on 
a substantial negative power of the temperature, the ther- 
mal conductivity improves for lower temperature, in the 
sense that a shallower gradient is needed to carry a given 
heat flux. Between T = 10^ K and 10^ K the radiative 
loss function /(T) ~ and this decrease in the efficiency 
of radiative loss at lower temperature, along with the im- 
provement in K± there, means that the temperature gradi- 
ent is very shallow at the cold outer edge of the shell, and 
steepens particularly sharply at the heated inner bound- 
ary, quite unlike the profile in a normal gas or solid. 



3 Analytical temperature profiles 

The function f{T) has an intricate dependence on temper- 
ature, but for our purposes we ma y take it to be a series 
of power laws (Rosner et al., 1978) 



In terms of the parameter 
2 



the heat flux density is 

Q{s) = koo^t' {s)/R. 



(8) 



(9) 



(10) 



With this form for Q{s), and the above form for /(T), we 
can rewrite Eqn. (|^) as 



— (st') -I- CST ^ = 0, 

as 



(11) 



in terms of a constant c = (3/2) ^x{RB)^ / ki, the param- 
eter b = 2{f3 — 2)/3, and boundary condition t'(1) = 0. 

This does not seem to be analytically solvable in general, 
but for the temperature range below 10^ K, we can set 
X = lO-'^'Wm^K-^ and /? = 2 in Eqn. (|). For /3 = 2, 
we have b — 0, and Eqn. (^ij) is easily solved to give 



:(1--^) 



(12) 



r = r{p) + ^(\n^^ + ^{p'-s^)], (13) 



3.1 Consequences of these solutions 

We may make a number of observations about these so- 
lutions. The first is that, if we define (5 = 1 — p, we can 
rewrite Eqn. ( p^ ) as 

{RSBni + OiS)) ^ ^(r(l) - T{p)), (14) 

and find that the radial thickness of the shell, RS, is, 
to 0{S), dependent only on the magnetic field and the tem- 
perature at the edge of the cell r(l) (where r(l) 3> t(p), 
or r(l) < Tp). For T(l) = lO^'K, Eqn. (|l|) pro- 
duces R6B « 0.12, giving a shell thickness of only 12 m 
for B = 10-2 T. 

The total power removed from a single core in length / 

is 

Pp = 2TTpRlQ{p). (15) 

If we have a large number of such filaments, each of ra- 
dius R, making up a total cross-sectional area of A, and 
we write Pa for the total power input across this area, 
then we can rewrite Eqn. (|l^) with the help of Eqn. ( |To| ) 
to find 

-^(1-A (16) 



lA 



This indicates two things. Firstly, since p has a mini- 
mum of zero, there is a maximum to the power that can 
be removed from the core by conduction into the shell 
(namely Pa = lO^^W for / = lO^m, A = W^xn^ and 
a = lO^^Km"'^) and this maximum is of similar order 
to the power of a large flare. It is difhcult to know if 
this is a coincidence or not: noting that it corresponds to 
the central core shrinking to a line, it seems likely that 
it merely reflects the observation that there is an upper 
limit to the power a given volume of plasma can radiate, 
and that Eqn. ( |l6|) is independent of n and T only because 
we have approximated /(T) oc to obtain the solution 
which produced it, making the losses n?' f{T) a constant. 
It is nonetheless true that n?f{T) is approximately con- 
stant below 10^ K, so that Eqn. (|l^) retains a qualitative 
general validity. 



with t{p) ^T{Tp). 



Secondly, Eqn. (16) is independent of the cell radius R. 
This is initially surprising, as one would expect that as 
the number of cells is increased and R falls, the ratio of 
the conducting surface area to the radiating volume would 
increase with 1/R and a greater power could be removed 
for a given value of p. However, since r' scales like R^, the 
flux density Q{p) itself scales like R, and so the total flux 
removed from the core scales with the same power of R as 
the radiating volume. 
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What is happening is that the temperature gradient at 
the surface of the core becomes steeper as the volume of 
shell it has to supply increases, and this steepening al- 
lows us to put limits on the permissible values of R. If 
the temperature gradient becomes too great, then normal 
conduction will be replaced by a more efficient anomalous 
mechanism which will effectively reimpose the maximum 
gradient allowed by normal conduction - the plasma will 
have this gradient for some distance from the surface of 
the core, until the solution Eqn. (|3|) takes over, with a 
suitably modified inner boundary condition, but the same 
outer one. The discussion below formally considers limits 
on the model's validity, but since the maximum-gradient 
solution is unlikely to be stable, it describes physical limits 
as well. 

To establish the limits on the radius in the absence of 
this effect, we must insist first of all that the shell width 
is larger than the ion-gyroradius, 



Figure 1 is in file figrl.ps 



R5 > Rr, 



where 



Rv 



10 



1/2 



B 



(17) 



(18) 



(in the numerical expressions here and below, B is in Tesla, 
R in metres, and n in m~^). Comparison with Eqn. ( |l^ 
shows that Eqn. (|l7| ) is comfortably satisfied, indepen- 
dently of B, for T < 10^ K. Secondly we can define a 
length-scale for the gradient. It, via dT/dr ^ Tp/lT so 
that, using the solution in Eqn. (|lj) and the definition 
of T in Eqn. (0), we have 



2 X lO'^p 
Tp^^RB^ ' 



and the condition It ^ Rp results in 



i?< 2 X 10^ 



BT^- 



(19) 



(20) 



Figure 1: The restrictions on the cell radius R, obtained 
by comparing the temperature gradient with the ion gyro- 
radius and ion mean free path. See text for explanation. 
For comparison, line 5 corresponds to the scale 1/fcc in 
Eqn. (Eq) below. 



Line 1 corresponds to Eqn. (pOD, line 2 to Eqn. (22), and 
lines 3 and 4 to the constraints R> Rp and R> Rf. All 
these are for the very weak field B = IQ-"^ T; the dashed 
lines show the corresponding constraints for B = 10~^ T. 

Finally, we may use the solution in Eqn. (|l^) to find the 
temperatures at s = 1 which are consistent with T{p) = 
Tp = 10^ K. These are displayed in Fig. ||. As can be seen 
there, constraining r(l) to lie between lO'' and 5 x 10'* K 
puts a tight constraint on the possible values of {p,RB). 
That is, only a narrow range of these parameters allow 
solutions with a cool outer shell in this T range. 



A third length scale is provided by the ion mean free path ^ ^ Anomalous enhancements to 



Rf 



2 X 10* 



(21) 



independently of the magnetic field. Imposing It Z Rf 
results in 

-,19 P 



i? < 10^ 



B^T, 



9/2 " 



(22) 



Together, these various constraints require that the ra- 
dius R lies between the solid lines (1) to (4) in Fig. |l|. 



The length scales in this model are consistently very small, 
due to the extreme smallness of the perpendicular con- 
duction coefficient ki in Eqn. (^. Observations in toka- 
maks consistently show thermal (and mass) transport co- 
efficients enhanced by one or two orders of magnitude over 
the classical values. This is generally attributed to micro- 
turbulence (small-scale, low-frequency, and more likely in 
the magnetic than the electric field), or to some other 
small-scale quasi-diffusive process, though theoretical ex- 
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Figure 2a and b are in files figtaua.ps and figtaub.ps 



Figure 2: The values of r(l) which are consistent with 
a boundary temperature T{p) = 10^ K, for various values 
of the parameters p and RB, taken from Eqn. ( p^ (/3 ^ 
2). The contours shown on the surface are at lO^K and 
5 X 10"* K. Fig. ||b shows the same two contours for (3 — 2 as 
well as the numerical solutions of the cases /3 = 1.9 (dotted 
line) and the polynomial fit (dashed line). Finally, the dot- 
dash line shows the contours for the case B = constant, 
discussed in Sect. H below. 



planations are inconclusive (Liewer, 1985). Though we 
should, of course, be wary of exporting these observations 
from the laboratory to the more tenuous and cooler at- 
mosphere of the sun, it seems safe to assume that the 
scales derived here are substantial underestimates. Ef- 
fective enhancement of ki will proportionately enhance 
the right-hand sides of Eqns. (pO|), ( p2[ ) and (|l4) and 
raise the corresponding upper boundaries in Fig. U In 
the Sun, ion-sound turbulence can result in a collision 



time much smaller than that given above dPriest, 1984| ), 
but this occurs only when the field has scales below 
~ 102i(B/T)(ne/m-3)-i(re/K)-i/2j^^^ comparable with 

the lower limits given in Fig. 

If we take a coronal magnetic field of 10^^ T and a 
boundary temperature of Tp = 10^ K, then Fig. |l] and 
the above two-order enhancement in the length scales give 
a filament scale of 10^ m, and hence require of the order 
of 10^ filaments to fill up a flare area of 10^'^ m^. 

3.3 Numerical solution of the tempera- 
ture profile 

For comparison with the above analytic results, and to 
confirm that our analytic solution is not somehow patho- 
logical, we have also solved Eqn. (^ numerically, with 
/(T) represented by power-laws with various values of /3, 



and by a polynomial fit ( Keddie, 197C| ) to the Cox- Tucker 
curve. The solutions are dependent on /3, but not unduly 
so, so that we can be confident that our results are quali- 
tatively insensitive to the precise form of the loss-function. 
See Fig. |. 



3.4 Comparison with the longitudinal 
temperature profile 

In section || above, we stated that we can ignore the ef- 
fects of longitudinal conduction, and we can now check 
that this is indeed true. We can avoid solving the full 2- 
d problem in the (s, z) plane, and instead calculate what 
the longitudinal temperature gradient would have to be 
for this assumption to fail. Writing the longitudinal flux 
density Q\\{z) as 



gii(z) = Kj|r5/2r, 



(23) 



where the 
IQ-^Wm-iR- 



coefficient 
^/2)/lnA dSpitzer, 



"11 



(2.23 

and T, 



dT/dz, we can use the solution in Eqn. (|12|) to write the 
ratio of transverse and longitudinal flux densities as 



Q± 
Q\\ 



2k° 



s 



-r{s)] — , 
2 ^ V T ' 



(24) 



where t(s) is given by Eqn. (p^). 

In Fig. ^, we show Tz{s), deflned so that Qj_ > Qn for 
Tz < Tz, for a variety of radii and magnetic field strengths. 
Note that Tz decreases towards the surface of the core 
s = p: the increase in temperature here enhances the lon- 
gitudinal conduction coefficient as it suppresses the trans- 
verse one. Note also that Tz increases with increasing R: 
this is initially surprising, as one would expect the as- 
sumption Q± ^ Q\\ to be more comfortably satisfied with 
finer filamentation; instead this reflects the steepening of 
the transverse gradient with larger R, which was noted 
above Eqn. (p^). In all cases, Tz{s = 1) = 0, reflecting 
the boundary condition on Eqn. ( [Til ) of dr/ds|s=i = 0, 
and we should also take Tz{s — p) = since, by hypoth- 
esis, the temperature at the surface of the core is a con- 
stant (Eqn. (|^)). It is difficult to estimate Tz in a model- 
independent way, but if we take Tz — (10^K)/Z ~ 1 K m~-^, 
we can see that for this geometry, the transverse flux den- 
sity should dominate over the longitudinal one everywhere 
except at the edge of the shell, s = 1. 
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part of our analysis, but the latter are more appropriate 
for perturbations with timescales shorter than the sound 
travel time. 



4.1 Isobaric perturbations 



Figure 1 is in file w7173.ps 



Figure 3: Values of Tz{s)/Km^^ , the values of the lon- 
gitudinal temperature gradient Tz{s) which make the 
transverse and longitudinal heat flux densities equal, as 
obtained from Eqn. (|2^). The curves are labelled by 
(logio(i?/m),logio(i?/T)). 

4 Stability 

We have considered only a static solution to the energy- 
balance problem. To determine the stability of these so- 
lutions, we must examine the heat equation 



1 dp J p dn 



7 — 1 dt 7 — 1 n dt 
The flux term K{n,T) is 
Kin,T) = 



+ K{n,T)+L{n,T) = 0. (25) 



1 d 

r ar 



dT 
dr 



and the radiative loss function L{n,T)^ ii?f{T), so that 
rK{n, T) + rL(n, T) = is our Eqn. (0). 

We combine this with the equation of state, the equation 
of motion, and the mass continuity condition, and perturb 
the equilibrium T, p and p by terms of the form 

a{r, t) — a^) + ai exp(i/t -f ikr). 

For simplicity, we will separately consider the limiting 
cases of isobaric [pi = 0) and isochoric {pi = 0) pertur- 
bations - the former match the conditions in the earlier 



For isobaric perturbations, pi = 0, the linearised heat 
equation, Eqn. (bsh, becomes 



7 



v—ni — g'^koTq ^^^i/cTi [1 /r + ik] 



7-1 no 



LtTi + LnUi = 0, (26) 



where Lt and _L„ denote partial derivatives of the heating 
function L{n,T). Combining this with the equation of 
state a = nT, the linearised perturbed equation of state 



rii 
no 



Ti 
To 



0, 



the pressure po — ksa (where fcs is Boltzmann's con- 
stant), setting /(T) = xT^ ^ and setting 7 = 5/3, gives 
the real part of the growth parameter v as 



Re V 



2 

5 ks 



a ( 



-3/2i,2 



(2 - (i)xT'o 



13-2 



(27) 



which is negative when (3 < 2 and k > k^ 
wavenumber fcj, defined by 



(2 



for a critical 



(28) 



Thus, solutions to Eqn. (^ are stable to isobaric perturba- 
tions with a scale shorter than l//cc- This applies also to 
more general forms of the loss function f(T), which can be 
locally approximated by a power-law . The dependence 
of 1/fcc on T and B is shown in Fig. |l]. 

4.2 Isochoric perturbations 

The isochoric case pi = is more appropriate for pertur- 
bations which are rapid compared to the sound time. The 
same analysis as above produces a growth rate 



1 



Po 



Rev ^ -nQ{f3xTo ^ + koT^ 



1- 



I To 



(29) 



which is certainly negative as long as (3 is positive. 

We can conclude that the only unstable perturbations 
are long-wavelength isobaric ones, which do not gener- 
ate sufficiently steep temperature gradients, though these 
perturbations should be stabilised when the perturbation 
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has heated enough to provide a sufficiently large gradient. 
Note that the solution outlined above Eqn. (pT|), where the 
temperature gradient in the inner part of the shell would 
have a constant maximum value, would be unstable also 
- since any local temperature enhancement will not result 
in a stabilising increase in the gradient, the temperature 
increase will go unchecked. 



5 Results in low-beta plasmas 



The full pressure-balance condition is 



n{s)T{s) + 



= a', (constant) 



(30) 



which we reduced to Eqn. (^) by taking the magnetic 
field to be constant across the filament. The two terms 
in Eqn. (|o|) are equal (for a = lO^^Km^^) when B sa 
1.9xlO~'^T. Ifi?is small compared with this, the condi- 
tion in Eqn. ( ^o| ) reduces to Eqn. (H). For a low-/? plasma, 
on the other hand, a change in the gas pressure can be 
compensated by a tiny change in the magnetic field, and 
Eqn. ( |30| ) reduces to i? = constant (we thank the referee. 
Dr. Schiifiler, for emphasising this point). 

In this latter case, given that the filament length and 
total flux are constant, and the field lines are frozen to the 
plasma, we can further take B/n = constant, and write 
Eqn. (0) as 



2/3 



(31) 



6 Conclusion 

We have presented a simple model of heat transport across 
the magnetic field in a highly filamented solar flare. This 
is analytically soluble in a simple case, which numerical 
solution shows to be representative. 

If we assume that the various potential complications 
do not qualitatively change the results, then the analysis 
very directly leads to predicted maximum cell sizes (Fig. |l|) 
which, though smaller than any yet observed, are larger 
than the lower limits imposed by the plasma conditions. 
Current uncertainties in our understanding of anomalous 
enhancements to the cross-field conduction strongly sug- 
gest that our upper limits on the cell size are underesti- 
mates. 

The upper limits we have derived here are strongly de- 
pendent on temperature and field strength, so that in the 
cooler and less magnetic environment of prominences, we 
would expect this mechanism to have a more prominent 
role than in fiares. 

Given that the mechanism described here is important 
in fiares of all sizes, then the dependence of p on Pa in 
Eqn. (^6|), and the interdependence of p, and RB in Fig. || 
are testable predictions. Even if we do not believe that this 
mechanism governs flare geometry as generally as this, we 
must regard the scale limits as delimiting a domain within 
which the effects of cross-field thermal conduction domi- 
nate over longitudinal conduction, and cannot be ignored. 
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